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Abstract 

A recent investigation has led to the possibihty of the existence of an interesting 
region of the asymptotic quasinormal mode spectrum of Schwarzschild-anti de Sitter 
black holes. In this asymptotic region, the real part of quasinormal mode frequen- 
cies approaches infinity while the damping rate approaches a finite value. These 
quasinormal mode frequencies were calculated using an analytic technique based on 
the complex coordinate WKB method. In this paper, we use a different analytic 
technique to calculate such quasinormal mode frequencies. The results of this paper 
provide further support of the possibility of the existence of these modes. 



The asymptotic quasinormal mode (QNM) spectrum of four dimensional Schwarzschild- 
anti de Sitter (AdS) black holes were calculated analytically for the first time by Cardoso 
et a/ in |1] , where the authors developed a method based on the monodromy technique of 
Motl and Neitzke[2l |3]. Later, Natario and Schiappafi] generalized the analytic results of 
|lj to spacetime dimensions greater than four. In a recent paper [5], Daghigh and Green 
reanalyzed the asymptotic QNM spectrum of Schwarzschild-AdS black holes in arbitrary 
spacetime dimensions [D > 4) using a different analytic technique called the complex co- 
ordinate WKB method O 0, E]- This method was first applied to asymptotic QNMs of 
black holes by Andersson and Howls in [9J, where the authors calculate the asymptotic 
QNM frequencies of non-rotating black holes. The authors of [5J not only confirm the 
analytic results found by Cardoso et al^ and Natario and Schiappa|4j, but also demon- 
strate that in certain spacetime dimensions the analytic techniques they use predict the 
existence of new regions of the asymptotic QNM frequency spectrum. These regions have 
not previously appeared or been explored in the literature [TOl [HI [121 [13 [IH [IS]- The 
QNMs of one of these regions, which appears only in even spacetime dimensions, resembles 
the normal modes in a pure AdS space. The real part of the QNM frequencies of such 
a region approaches infinity while the damping rate approaches a finite value. In [S], the 
name "highly real" is chosen for these modes to correspond to the widely used term "highly 
damped" , which is the situation when the imaginary part of the frequency is much larger 
than the real part. If the highly real QNMs do indeed exist, it is argued in |5] that these 
modes become the most relevant modes in the context of AdS-CFT correspondence. They 
can also serve as a testing ground for the new interpretation of the black hole asymptotic 
QNMs proposed by Maggiore|,16j. The very interesting consequences of the existence of 
highly real QNMs motivates numerical and analytical investigations to prove or disprove 
the existence of these modes. 

In this paper, we use the analytic technique which was developed by Cardoso et al in 
[1] to see if we can reproduce the results of [5] for highly real QNMs. At first glance, it 
is not obvious whether the two different analytic techniques mentioned above produce the 
same results. Using the complex coordinate WKB method, Daghigh and Green [5] were 
able to obtain two unphysical and two physical solutions for highly real QNM frequencies 
by changing the dominancy of the WKB solutions on the pairs of Stokes lines, which 
are crossed on the path taken along anti-Stokes lines, and by rotating clockwise (in the 
upper half of the complex plane) and counter-clockwise (in the lower half of the complex 
plane) close to the singularity at the origin of the complex plane (r = 0). In the method 
developed by Cardoso et a/[lj, there are no use of Stokefl lines. Therefore, one can only 
find two solutions for the highly real QNMs by rotating clockwise (path 1 in Fig. [1]) or 
counterclockwise (path 2 in Fig. [1]) close to the origin of the complex plane. These two 
solutions could be the unphysical ones which were discarded by the authors of [5]. If that 
is the case, then we can conclude that we have found strong evidence against the existence 
of the highly real QNMs. In what follows, we derive an expression for highly real QNM 
frequencies using the analytic technique developed in pQ. 

Various classes of non-rotating black hole metric perturbations in a spacetime with 

^Note that what we call anti-Stokes lines in this paper are called Stokes lines in [Ij and [4]. 
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dimension D > A are governed generically by a Schrodinger wave-like equation of the form 

^ + K-^[rW]}^ = 0, (1) 

where V{r) is the QNM potential obtained by Ishibashi and Kodama[T71 [181 [IH] for scalar 
(reducing to polar at D = 4), vector (reducing to axial at D = 4), and tensor (non- 
existing at D = 4) perturbations. It is worthwhile to mention that the effective potential 
for tensor perturbations is equivalent to that of the decay of a test scalar field in a black 
hole background in every spacetime dimensions including D = A. Following [T] and [4], we 
assume the perturbations depend on time as e*"^*. Consequently, in order to have damping, 
the imaginary part of a; (a;/)must be positive. Keep in mind that in |^J, the authors assume 
the perturbations depend on time as e"""^* instead of e*"^*. In Eq. ([1]), the tortoise coordinate 
z is defined by 

dr , , 

where /(r) is related to the spacetime geometry, and is given by 

m^'-^-^-r'- (3) 
The ADM mass, M, of the black hole is related to the parameter /i by 

where Gd is the Newton gravitational constant in spacetime dimension D and is the 
area of a unit n-sphere, 

o n+i 
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The value of the cosmological constant. A, is given by 

The effective potential in AdS space is zero at the event horizon (z — oo) and diverges 
at spatial infinity [z — > 77, where is a constant discussed in detail in [T| 15)^1 ) . Therefore, 
the asymptotic behavior of the solutions is 

T / \ I 6 cLS Z ^ OC / ^\ 

which represents an "out-going" wave at the event horizon and no waves at spatial infinity. 

To determine the condition on highly real QNMs using the method developed by Car- 
doso et a/[lj, we need to know the behavior of anti-Stokes lines in the complex r-plane. 



^The symbol used in [I] and [J for 77 is xq. 
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Figure 1: Schematic presentation of the anti-Stokes hnes in the complex r-plane for 
Schwarzschild-AdS black holes in four spacetime dimensions. Rh is the event horizon. 
The two dashed arrows show the paths we take along anti-Stokes lines to determine the 
condition on the highly real QNM frequency. 



The details of this process have been explained in detail in [H H]. As an example, in 
Fig. [U we show the schematic behavior of the anti-Stokes lines in the complex r-plane 
for Schwarzschild-AdS black holes in D = 4, when 77 is a negative real constant. For large 
black holes ^ 

where Th is the Hawking temperature. In spacetime dimensions six and higher, we will still 
have an anti-Stokes line that extends to infinity along the negative real axis and an anti- 
Stokes line that extends to the event horizon, Rh, along the positive real axis. The only 
difference is that in higher dimensions we have more anti-Stokes lines extending above and 
below the real axis. In order to determine the condition on highly real QNM frequencies, the 
same paths shown in Fig. [1] will be taken in every spacetime dimension. For numerically 
generated pictures of anti-Stokes lines in D > 4, refer to [5]. Once we determine the 
structure of anti-Stokes lines, we need to know the behavior of the wavefunction \1/ near 
the singularity at r = and when r — 00. Near r = 0, the tortoise coordinate simplifies 
to 



j.D~2 



2(D-2)/i 



(9) 



and the potential simplifies to 



y vi^)] - 4;^ ' (10) 

where j = 2 for vector perturbations and j = for scalar and tensor perturbations [1]. In 
this region, the differential equation ([1]) can be solved exactly. The solution is 

-^{z) = A+V2ttuuzJj_ (ujz) + A_V27wjzJ_2. (uz) , (11) 
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where Jy represents a Bessel function of the first kind and A± are complex constants. 
Note that D is an even number because highly real QNMs only exist in even spacetime 
dimensions [5]. Therefore, it is clear that near the origin of the complex plane, ujz <^ — 1 
on the anti-Stokes line that extends to —oo along the real axis. As a result, we can use the 
approximation 

/ 2 Z/TT TT 

Ju[x) ^ \ — cosfa; H 1 — ) when x ^ — 1 (12) 

V TTX 2 4 

to find that 

^(z) ~ 2A+cos{ujz + a+) + 2A^cos{loz + a^) 

= {A+e'"+ + A_e^°-)e^'^^ + {A+e-'"+ + A_e-'"-)e-'''' . (13) 

Here 

a± = ^(l±j). (14) 

The result (1130 is valid on any anti-Stokes line with negative uz. Note that in addition to 
the anti-Stokes line on the negative real axis, uz is negative on the anti-Stokes line along 
the positive real axis. 

In the limit where r oo, the tortoise coordinate simplifies to 

1 , , 

Z^T]-—-. 15 

|A|r 

In this limit, the QNM potential also takes the simple form 

where joo = D — 1,D — 3,D — 5 for tensor, vector, and scalar perturbations respectively. 
Note that the potential (fT6|) is zero for vector (axial) and scalar (polar) perturbations 
in four spacetime dimensions. In reality, the QNM potential for these perturbations in 
four spacetime dimensions asymptotes to a constant as r — oo (z — ^ ?^). This may 
look problematic, but note that we are doing the calculations in the asymptotic limit 
where |a;| oo. Therefore, in the wave equation ([T]), replacing the QNM potential, which 
asymptotes to a constant for large r, with a zero should not in principal affect the outcome. 
In fact, Natario and Schiappa in [4J have used exactly the same potential in Eq. f|T6l) and 
their results for polar and axial perturbations in four spacetime dimensions are in perfect 
agreement with the previous numerical calculations. It is clear that such modification of 
the potential at large r does not affect the gap term in the QNM frequency. However, 
we suspect that this modification may affect the offset term. In fact, it is indicated in [1] 
that the offset term of polar perturbations does not match the numerical results in [TT] for 
large black holes. This mismatch may originate from the modification of the potential in 
four spacetime dimensions. This issue clearly needs further investigation. The differential 
equation ([1]) can be solved exactly in the region r ~ oo and the solution is 



^{z) ~ B+J2tiuj{z - ri)Jj^ {uj{z - r])) + B_J2nu{z - r])J-i^ {^{z - v)) , (17) 



2 V ^ - 2 
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where B± are complex constants. The boundary condition at infinity tells us that the wave 
function \1/ should die off as r oo. This requires 5_ = 0. Also, since 

uj{z-v)---7^ (18) 
|A|r 

has a positive real value on the anti-Stokes line along the negative real axis, we can use 
the approximation 

/ 2 lyiT 7T 

Jyix) ~ \\ — cos(x ) when x ^ 1 (19) 

\ Tlx 2 4 

to write the solution (JTTll as 

■^{z) ~ 25+ cos[u;(z -r])~l3]= S+e-^^e^'^"-") + S+e^^e-^'^^"-^") , (20) 

where 

/?=^(l+Joo). (21) 

Note that the wavefunction \1/ should not change in character along anti-Stokes lines. 
Therefore combining Eqs. (IT^ and (1201) results in the first condition on ^4+ and A_, where 

(^+e^°+ + A_e^°-) e^'^e^^ = (^+6"^°+ + A.e"*"-) e-^^^e"*'^ . (22) 

In order to impose the boundary condition at the event horizon, we need to move to the 
anti-Stokes line on the positive real axis which ends at the horizon. We can do this by a 
rotation of either — tt (path 1 in Fig. [1]) or +7r (path 2 in Fig. [1]) near the singularity at 
r = 0. The rotation of ±7r in the r coordinate corresponds to a rotation of ±7r(L' — 1) in 
the tortoise coordinate. From 

J^{x) = xy{x) , (23) 

where (p is an even holomorphic function of x, it can be shown that after a rotation of 
—n{D — 1) in the tortoise coordinate one can write 

A±v^^±2(^-2) ~ e"'^(^"^)°± cos(a;;z + a±) when CJ2 < -1 . (24) 

Therefore, after a rotation of —7i{D — 1), the solution (1131) changes to 

^(2) ~ 2A+e-*2(^5-2)a+ pQg^^^ ^ ^ 2A_e-*2P-2)"- cos{ujz + a_) 

+ (A+e-'[2{D-2)+l]a+ ^^_g-^[2{D-2)+l]a_j g-i.;. _ ^25) 

According to the boundary condition (JTj), the solution on the anti-Stokes line ending at 
the event horizon should behave as e*"^^. This fact results in a second condition on the 
coefficients and A^, where 

^^g-^[2{D-2)+l]a+ ^ ^_g-i[2{D-2)+l]a_ _ Q _ (26) 
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Equations fl22|) and fl26|) can only have nontrivial solutions if 

g-i[2(D-2)+l]a+ g-i[2(D-2)+l]a_ 



-ia+ 



-iLori 



g i/3 gia_ 



-iujrj p—if) 



(27) 



In four spacetime dimensions, this equation leads to 



g2iaj77g2i/3 



g— i5Q!+ g— i«- 



Q—iba+ Qia— g— i5a_gjQ+ 

sin[3(a+ — a_)] 
2cos(a+-a_ 



4 cos^(a+ 



1 



(2^ 



Since + a_ = | and q;+ 



|7r, it is easy to show that the above equation leads to 



TT 2 

cur? = —B ~ niT 1 — In 

' ^ 4 2 



2 cos -TT 

2 



1 



2 cos 



as n 



oo 



(29) 



This result is in perfect agreement with the result obtained by Daghigh and Green in 
[S]. Of course, the sign difference in the imaginary part of wq is due to the fact that our 
perturbations depend on time as e*"^* while Daghigh and Green use e~*^*. Since 77 is a 
negative real constant according to both pLj and Eq. (129!) leads to QNM frequencies 
with negative imaginary part for all types of perturbations. QNMs with uj < cannot 
be physical because these modes grow with time instead of getting damped. By repeating 
the same calculations starting from Eq. (!27|l . it is easy to show that, in D > 6 spacetime 
dimensions, we also obtain QNM frequencies with negative damping which cannot be 
physical. 

Let us now consider path 2 in Fig. (H where we have to make a rotation of vr instead 
of — vr in the r coordinate. We can repeat the steps through (I27p with a rotation of 
7r(D — 1) in the tortoise coordinate to find 



J[2{D-2)-l]a+ 



-ia+ 



i[2(D-2)-l]o_ 



-if3 



. 



(30) 



This equation leads to a condition on u, where 



g2«i^-'?g2i/3 



i[2{D-2)-l]a+^ 



J[2{D-2)-l]a- 



gi[2{D-2)-l]a+gja_ _ gj[2{D-2)-l]a_ gj 

.-»fa^+a-) Sin[(D - 2){a+ - a_)] 
sin[(i5 — 3)(«_|_ — a_)] 



-i{a++a^)£ — ^cos(a+ -a_ 



D-3 

-22 . 2 
— sm 



(a+ - «_) + \i^-^f-niD-2f-^^] gi^4^ 



sm [a^ — a^) — 



1 - 



3! 



-1^ • 2/ 

— sm (a+ 



(31) 
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After entering the values for a± , it is easy to show that the above condition on u gives 



This resuh again matches perfectly with the result obtained by Daghigh and Green in [5] 
considering the difference in the choice of the time dependency of the perturbations. Since 
?7 is a negative real constant, it is clear that these QNM frequencies have positive imag- 
inary part. This indicates the presence of damping which makes these QNMs physically 
meaningful. For more arguments on why these modes may be physical see [5]. 

Conclusions 

In this paper, we use a different analytic technique to confirm the results obtained by 
Daghigh and Green in [S]. The two analytic methods lead to identical solutions. The 
only difference is that the complex coordinate WKB method which is used in [5J gives 
two physical and two unphysical solutions, but the technique used in this paper results 
in only one physical (Eq. fl32|) ) and one unphysical (Eq. fl29l) ) solution. No difference has 
been observed between the outcome of these two analytic methods in the past. Further 
investigation of the difference found in this paper may lead to a deeper understanding 
of these two analytic techniques and their possible limitations. While this difference is 
interesting, we do not see it as a major concern. After all, the calculations are done in 
the asymptotic limit where the overtone number, n, approaches infinity. In the large n 
limit, we can neglect all the finite terms in the real part of the highly real QNM frequency. 
Neglecting the finite terms will resolve the difference between the two analytic techniques. 
The results of this paper leave no doubt that the analytic techniques in the literature [H [5] 
predict the existence of the highly real QNMs of Schwarzschild-AdS black holes. If these 
modes do not exist, the big question is why they appear in analytic calculations which are 
very successful in producing the correct asymptotic QNM frequencies of Schwarzschild-AdS 
black holes. If these modes do indeed exist, then the question is why the previous numerical 
calculations[ini[IIl[I21IIl] were not able to detect such modes. It is important to note that 
numerical calculations should be able to generate the unphysical QNM frequencies as well 
as the physical ones (see for example [20]). In order to shed light on these issues, further 
numerical and analytical investigations are necessary. 

Acknowledgments 

We are grateful to Michael Green for numerous useful discussions on the details of these 
calculations. 



References 

[1] V. Cardoso, J. Natario, and R. Schiappa, J. Math. Phys. 45 (2004) 4698. 

[2] L. Motl, Adv. Theoret. Math. Phys. 6 (2003) 1135. 

[3] L. Motl and A. Neitzke, Adv. Theoret. Math. Phys. 7 (2003) 307. 

[4] J. Natario and R. Schiappa, Adv. Theoret. Math. Phys. 8 (2004) 1001. 




as n 



oo. 



(32) 



8 



[5] R.G. Daghigh and M.D. Green, larXiv:0808. 1596 [gr-qc]. 

[6] J. Heading, An introduction to phase-integral methods (Methuen, London, 1962) 

[7] N. Froman and P.O. Froman, JWKB approximation: contributions to the theory (North- 
Hohand, Amsterdam, 1965) 

[8] M.V. Berry and K.E. Mount, Reps. Prog. Phys. 35 (1972) 315. 

[9] N. Andersson and C. J. Howls, Class. Quant. Grav. 21 (2004) 1623. 

[10] E. Berti and K.D. Kokkotas, Phys. Rev. D67 (2003) 064020. 

[11] V. Cardoso, R. Konoplya, J.P.S. Lemos, Phys. Rev. D68 (2003) 044024. 

[ 12 ] R.A. Konoplya, Phys. Rev. D66 (2002) 044009. 

[13] S. Musiri, S. Ness, and G. Siopsis, Phys. Rev. D73 (2006) 064001. 

[14] A.S. Miranda and V.T. Zanchin, Phys. Rev. D73 (2006) 064034. 

[15] A. Ghosh, S. Shankaranarayanan, and S. Das, Class. Quant. Crav. 23 (2006) 1851. 

[16] M. Maggiore, Phys. Rev. Lett. 100 (2008) 141301. 

[17] A. Ishibashi and H. Kodama, Prog. Theoret. Phys. 110 (2003) 701. 

[18] A. Ishibashi and H. Kodama, Prog. Theoret. Phys. 110 (2003) 901. 

[19] A. Ishibashi and H. Kodama, Prog. Theoret. Phys. Ill (2004) 29. 

[20] E.W. Leaver, Proc. Roy. Soc. Lond. A 402 (1985) 285. 



9 



